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compared. The current capability for these two instruments to measure and discriminate ICI 
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I. Introduction 


HE ingestion of ice crystals into jet engines has been attributed to uncommanded jet engine power loss events 
during flight'>, and has become a significant focus for research in the atmospheric icing community. As a result, 
NASA has been conducting ice crystal icing (ICI) studies in the the NASA Propulsion Systems Laboratory (PSL) to 
advance the community’s understanding of the fundamental physics behind this aviation safety hazard®’. The unique 
nature of ICI conditions has slowed cloud characterization in PSL‘, but it has also stimulated the development of 
newer instrumentation to measure such conditions. The Artium Technologies, Inc. Phase Doppler Interferometer (PDI) 
and High Speed Imaging (HSI) instruments are two such examples that have been further developed as a result of the 
community’s interest in ICI and funding support through the Small Business Innovative Research (SBIR) program. 
This paper presents select cases of particle measurements acquired using the PDI and HSI instruments in PSL 
during the first Fundamentals of Ice Crystal Icing Physics Test conducted in March of 2016. The results presented 
herein examine the data acquired from the PDI and HSI during the test, comparing key parameters and evaluating the 
current capabilities of both instruments. The objective of this work is to increase understanding of the capabilities of 
these two instruments thereby allowing for improved characterization of mixed phase and ice crystal icing conditions. 


II. Experimental Description 


The first Fundamentals test generated mixed-phase ice accretions on a test article, a NACA 0012 airfoil wing 
section with a 267 mm chord, in a characterized cloud. The goal of the effort was to evaluate if PSL could be used for 
further Fundamentals tests by simulating conditions considered similar to those expected within a jet engine core 
during ICI conditions. Several instruments were used to characterize the mixed-phase cloud including systems from 
Droplet Measurement Technologies (DMT), Science Engineering Associates (SEA), the National Research Council 
of Canada (NRC), and NASA in addition to those from Artium Technologies, Inc. Besides the PDI and HSI from 
Artium Technologies, Inc., the DMT Cloud Droplet Probe (CDP) and Cloud Imaging Probe - Grey Scale (CIP-GS) 
were the only other particle sizing instruments used during the test. The results from the DMT instruments will be 
covered in a separate conference paper, Ref. 7. 


A. Phase Doppler Interferometer and High Speed Imaging Instruments 

The Artium Technologies, Inc. PDI is a single particle counter using a flux sampling technique. The physical 
principles underlying the PDI have been well documented in numerous publications, including Ref. 9. The Artium 
PDI system splits a laser beam and focuses the two resulting coherent beams to a common point in space, creating an 
interrogation volume and generating a local interference fringe pattern. Particles passing through this volume will 
scatter the light, creating a Doppler burst signal as they pass the interference fringe pattern. The Artium Technologies, 
Inc. PDI measures this Doppler burst with three detectors at separate spatial locations. The resulting phase shift of the 
Doppler burst signals allows measurement of the spacing of the interference fringe pattern, which is used to determine 
particle size. Using three detectors provides a means to avoid phase ambiguity when the phase cycles past 360°, and 
also provides redundant measurements that are used to validate the signals, offering an indication as to whether the 
particles are quasi-spherical or irregular-shaped ice particles. The system can size spherical and quasi-spherical 
particles. Irregularly shaped and partially melted particles will not follow the calibrated phase shift-size relationship, 
which provides an avenue for particle material phase and morphology discrimination. 

The Artium Technologies, Inc. HSI is another particle sizing instrument that uses a spatial sampling technique. 
The HSI acquires high-resolution images of particles passing through the interrogation volume. This volume is created 
by focusing several laser beams on a common spatial point, illuminating particles for image capture by a CMOS 
camera with a long range Schmidt—Cassegrain telescope, which is recording at 300 Hz. The lasers are simultaneously 
pulsed with a pulse duration on the nanosecond time-scale, reducing exposure times and motion blur of the images. 
With knowledge of the system resolution determined during setup calibration, the system can size spherical and 
irregularly shaped particles, and quantitative assessments of particle morphology can be made, allowing for potential 
identification of glaciated particles. The resolution for the HSI during the test was 1.8 m/pixel. 

The limitations of any drop sizing instrument should be understood when comparing systems. Thus, it should be 
noted that the HSI has a larger minimum measurable particle size than the PDI, which can measure particles below | 
tum, because a single darkened pixel is insufficient to identify and measure a particle in an image. Additionally, the 
size of the sample volume for the HSI is, in general, smaller than the PDI, and decreases in relative size to the PDI 
sample volume with increasing airspeed. This is because the sample volume size for the HSI is independent of the 
airspeed, while the sample volume size for the PDI is proportional to the airspeed. 
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The PDI and HSI used during the test were modular, and the non-intrusive nature of the instrumentation allowed 
the systems to be used across several test days without interfering with other measurements and studies. Figure | 
shows the experimental setup for the NACA 0012 ice accretion studies, which was used for 2 of the 8 testing days. 
The PDI and HSI, which are labeled in Fig. 1, were focused approximately 152 mm above the centerline for these 
days. The PDI and HSI were also used during the test day dedicated to the DMT CDP and CIP-GS particle 
measurements, but only during the CDP portion of testing. The PDI and HSI were focused approximately 152 mm 
above the centerline for this portion of the test as well. More information on the positioning of the instrumentation 
and the test article during the test is available in Ref. 7. Unfortunately, no measurements were made at the centerline 
of the duct with the PDI or HSI during the first Fundamentals test. 


B. Select Cases 

Table 1 outlines the select cases examined in this paper. The cases are separated into two series. The first series is 
a humidity sweep of the PSL plenum conditions at a large dyo.50, 50 um, based on the NASA Icing Research Tunnel 
(IRT) drop size calibration, and the second series is a humidity sweep at a small d,o.50, 15 um, also based on the IRT 
drop size calibration. Varying the plenum humidity varies the local wet bulb temperature, Twz. The wet bulb 
temperature affects the freeze-out of the cloud particles in the plenum, which is where the cloud is expected to have 
the greatest residence time’®. The test points in Table | are arranged in order of decreasing Twp. 

The station 1 total water content, TWCs,, values defined in Table | is are bulk averaged values, which assume that 
the cloud is uniformly distributed within the innermost 610 mm diameter of the station 1 plane. The spray bar nozzle 
pattern was optimized to generate a cloud that did not extend outside the innermost 610 mm diameter at station 1. 
Refer to Ref. 7 for more detail on the station 1 total water content and cloud uniformity. 


Table 1. Select Test Points 


Test Point PPL @PL Twe Uusi Ms; TWCs) dyo.s0t 
--- psi % °C m/s --- g/m} Mm 
Large dyo.so0 Series 
210 6.5 60 2.4 85 0.25 2 50 
209 6.5 47 0.6 85 0.25 2 50 
206 6.5 40 -0.4 85 0.25 2 50 
204 6.5 25 -2.6 85 0.25 2 50 
202 6.5 20 -3.3 85 0.25 2 50 
201 6.5 10 -4.9 85 0.25 2 50 
201.1 6.5 0 -6.0 85 0.25 2 50 
Small dyo.so Series 
211 6.5 60 2.4 85 0.25 2 15 
208 6.5 47 0.6 85 0.25 2 15 
207 6.5 40 -0.4 85 0.25 2 15 
207.1 6.5 38 -0.5 85 0.25 2 15 
205 6.5 25 -2.6 85 0.25 2 15 
203 6.5 20 -3.3 85 0.25 2 15 


* Values based on conditions in the plenum 
* Values based on NASA Icing Research Tunnel cloud drop size calibration and do not represent 
PDI or HSI test data 


III. Results and Discussion 


Figures referenced in this section are located at the end of the paper. The figures are arranged in the order of 
decreasing Twg as shown in Table | for each series, and each test point figure contains five subfigures, including (a) 
cumulative percent volume distribution comparison, (b) number density comparison, (c) water content comparison, 
(d) area-perimeter correlation, and (e) raw signal phase difference plots. While the first three plots provide useful 
comparison between the PDI and HSI, the area-perimeter and raw signal phase difference plots provide an indication 
of the morphology of the particles, which is useful when analyzing the results. The area-perimeter plots are developed 
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from the HSI data where spherical particles are expected to closely follow the theoretical line for a circle, AR = | line. 
The AR = 2 and AR = 3 represent ellipses with major diameters 2 and 3 times the unitary minor diameter, respectively, 
and are intended to provide frames of reference. Similarly, the raw signal phase difference plots are developed from 
the PDI data where spherical particles are expected to closely follow the dashed lines, as described in the Experimental 
Description section. 


A. Large dyo.so Series 

The PDI and HSI particle size distributions (PSD) agree reasonably well throughout the Large dyo,50 Series, which 
can be seen in Fig. 2 through 8. Examining Table 2, which shows normalized cumulative percent volume distribution 
values, the PSDs for both instruments appear independent of the decreasing Twz in the series. The number densities 
and water content distributions for both instruments also agree throughout the series, but appear low, which may 
indicate that the PDI and HSI were measuring close to the outside perimeter of the cloud. It is clear from the number 
density and water content comparisons that the HSI rolls off towards the lower end of the particle spectrum faster than 
the PDI. However, this behavior is expected because the HSI has a larger minimum measurable particle size than the 
PDI. At the other end of the spectrum, the HSI cuts off sooner that the PDI. Again, this behavior is expected because 
of difference in the sample volume size between two instruments. Throughout the rest of the distribution, the PDI and 
HSI agree well. 

Progressively examining each test point in the series reveals a gradual change in the raw signal phase difference 
plots from the PDI and the particle area-perimeter correlation from the HSI. As the Twz decreases, the scatter in the 
raw signal phase difference plots increases, and the data trends away from the theoretical line for a circle (AR = /) in 
the area-perimeter correlation, both indicating increasing presense of irregularly shaped ice particles. The sample 
particle images in Fig. 9 and 10 show the morphology difference between test points 210 and 201.1 supporting this 
indication. These items demonstrate the current ability to assess particle morphology using the PDI and HSI, and a 
future potential avenue for statistical morphology characterization of ICI conditions. The particles lying further from 
theoretical lines in both the raw signal phase plots and the area-perimeter correlation plots are believed to be glaciated, 
but further validation is required to corroborate such assessments. 


Table 2. Large dvo.so Series Normalized Cumulative Volume Distributions 


Test Point Amin dyo.10 dyo.25 dyo.50 dy0.90 dy0.99 max 
fhe ym pm pm pm pm Mm pm 
PDI 
210 6 26 37 53 97 137 164 
209 7 27 37 53 97 138 156 
206 5 27 37 53 96 139 155 
204 5 27 37 52 97 143 177 
202 4 27 37 53 97 130 141 
201 1 25 36 52 92 141 149 
201.1 1 26 36 52 94 139 143 
ASI 
210 6 24 33 47 98 129 137 
209 7 25 32 47 88 124 125 
206 5 25 33 47 89 159 160 
204 5 26 35 49 86 119 125 
202 4 26 36 50 90 116 130 
201 1 27 38 57 99 116 118 
201.1 1 28 37 ‘52: 87 117 123 


B. Small d,o.50 Series 

Inspecting Fig. 11 through 16, it is apparent the PDI and HSI PSD do not agree well in the Small do,50 Series. 
Acquiring the PSD with the HSI is difficult because the distribution is concentrated towards the lower end of the size 
spectrum, and the roll off of the HSI towards the lower end is more apparent in the Small d,o.50 Series number density 
and water content comparisons than the Large d,o.59 Series comparisons. Table 3 shows decreasing independence of 
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the PSD for the PDI to decreasing Twz, while the HSI PSD remains consistent, but shifted towards the higher end of 
size spectrum, throughout the series. This indicates that the PSD is not changing, but rather, the particles are becoming 
more irregularly shaped, reducing validation on the lower end of the particle spectrum for the PDI. These items 
indicate that further development of the PDI and HSI towards measuring the lower end of the irregular-shape particle 
spectrum is necessary. 

Again, each test point in the series reveals a gradual increase in the scatter in the raw signal phase difference plots 
from the PDI, similar to that observed in the Large dyo.50 Series. The data trending away from the AR = / line in the 
area-perimeter correlation plots from the HSI is evident upon inspection, but a bit easier to distinguish than the Large 
dyo.so Series. Finally, the sample particle images in Fig. 17 and 18 show the morphology difference between test points 
211 and 203, respectively. Thus, similar to the Large d,o.50 Series, the Small d,o.50 Series particle images show a clear 
increase in irregular-shape particle morphology with decreasing Twa, indicating the presence of particles believed to 
be glaciated. Again, further validation is required to corroborate such assessments. 


Table 3. Small d,o.s0 Series Normalized Cumulative Volume Distributions 


Test Point Amin dyo.10 dyo.25 dyo.50 dyo.90 dy0.99 Amax 
= ym pm pm pm ym ym ym 
PDI 
211 0.7 12 15 20 33 44 56 
208 0.6 11 15 20 33 47 56 
207 0.6 11 15 20 33 48 58 
207.1 0.5 12 16 22 43 65 82 
205 0.5 13 18 25 51 84 89 
203 0.5 13 18 26 54 83 87 
ASI 
211 16 21 25 32 42 45 46 
208 12 20 23 29 42 55 62 
207 11 19 23 28 41 55 62 
207.1 11 19 23 28 40 50 62 
205 11 19 23 28 41 50 55 
203 9 19 23 29 41 67 68 


IV. Conclusions 


NASA has completed the first Fundamentals of Ice Crystal Icing Physics test in PSL. The Artium Technologies, 
Inc. PDI and HSI instruments were successfully used to measure the icing cloud during the test. The PDI and HSI 
demonstrated good agreement during the Large d,o.s0 Series through a range of Twz. Due to minimum measurable 
range of the HSI, there was not good agreement between the PDI and HSI during the Small dyo.50 Series. Thus, further 
development of the PDI and HSI towards measuring the lower end of the irregular-shape particle spectrum is 
necessary. The PDI raw signal phase difference and HSI area-perimeter correlation plots demonstrated the current 
capability for the instruments to allow researchers to examine particle morphology, and consequently discriminate 
glaciated conditions. Further validation is required, but initial analyses appear to indicate a potential avenue for a 
future capability to statisitically characterize ice crystal cloud material phase with the PDI and HSI. Thus, the results 
from the test offer insight into the current capabilities of the PDI and HSI, and a representative data set of ICI 
conditions to drive future development towards improved measurement and characterization of hazardous icing 
conditions. 
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Figure 1. First Fundamental Ice Crystal Icing Test Configuration 2 setup, showing the Artium Technologies, 
Inc. Phase Doppler Interferometer and High Speed Imaging instruments focused approximately 152 


millimeters above the duct centerline 
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Figure 2. Test Point 210 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (e) PDI Raw Signal Phase Difference 
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Figure 3. Test Point 209 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (e) PDI Raw Signal Phase Difference 
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Figure 4. Test Point 206 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 


(d) HSI Particle Area-Perimeter Correlation , and (e) PDI Raw Signal Phase Difference 
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Figure 5. Test Point 204 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (e) PDI Raw Signal Phase Difference 
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Figure 6. Test Point 202 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (e) PDI Raw Signal Phase Difference 
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Figure 7. Test Point 201 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (e) PDI Raw Signal Phase Difference 
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Figure 8. Test Point 201.1 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (d) PDI Raw Signal Phase Difference 
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Figure 9. Test Point 210 sample particle images from the HSI 
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Figure 10. Test Point 201.1 sample particle images from the HSI 
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Figure 11. Test Point 211 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (d) PDI Raw Signal Phase Difference 
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Figure 12. Test Point 208 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (d) PDI Raw Signal Phase Difference 
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Figure 13. Test Point 207 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (d) PDI Raw Signal Phase Difference 
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Figure 14. Test Point 207.1 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (d) PDI Raw Signal Phase Difference 
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Figure 15. Test Point 205 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (d) PDI Raw Signal Phase Difference 
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Figure 16. Test Point 203 Results showing (a) Cumulative Percent Volume Distribution PDI-HSI Comparison, 
(b) Number Density Distribution PDI-HSI Comparison, (c) Water Content Distribution PDI-HSI Comparison, 
(d) HSI Particle Area-Perimeter Correlation , and (d) PDI Raw Signal Phase Difference 
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Figure 17. Test Point 211 sample particle images from the HSI 
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Figure 18. Test Point 203 sample particle images from the HSI 
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